Dynamical Coupling Between Ferromagnets Due to Spin Transfer Torque 
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We use a combination of analytic calculations and numerical simulations to demonstrate that 
electrical current flowing through a magnetic bilayer induces dynamical coupling between the layers. 
The coupling originates from the dependence of the spin transfer torque exerted on the layers on the 
relative orientations of their magnetic moments. We demonstrate that such coupling modifies the 
behaviors of both layers, significantly affecting the the stability of the current-induced dynamical 
regimes and the efficiency of current-induced magnetic reversal. 

PACS numbers: 85.75.-d, 75.60.Jk, 75.70.Cn 



Current-induced spin transfer (ST) effect |l| is the 
most promising mechanism for manipulation of magnetic 
nanodevices, due to the simplicity of the implementation 
and potential power benefits. The main obstacle for vi- 
able applications of the effect is the large magnitude of 
the required current, which is too close to the limit of the 
physical stability of devices. A basic magnetoelectronic 
device consists of a magnetic bilayer F1/N/F2, where Fi 
is a magnet needed to polarize the current, N is a metal- 
lic or insulating nonmagnetic spacer, and F2 is a nano- 
magnet whose magnetic configuration can be changed by 
current via ST. The efficiency of spin transfer can be 
characterized by the zero-temperature threshold current 
I c for the onset of current-induced magnetic dynamics of 
F2. In the framework of the widely accepted spin transfer 
torque (STT) model HQ, 



eH ef fm 2 a/(hg2) 



(1) 



where e is the electron charge, i? e // is the effective mag- 
netic field which includes the magnetic anisotropy of F2, 
a is the Gilbert damping parameter, and g 2 is a func- 
tion of the relative orientations of the magnetic moments 
mi and m 2 of Fi and F2, respectively, which depends 
predominantly on the spin polarizing properties of Fi. 

Several directions are pursued for reducing I c . The 
function g 2 is proportional to the polarization p = (If — 
of the current generated by Fi |4j. Here, If and 
If are the spin-up and spin-down contributions to the 
current. Therefore, I c can be reduced by enhancing the 
spin-polarizing properties of Fi. However, the difference 
between the typical value p 0.7 for the common Fi 
such as Py=Ni 80 Fe2o and the largest possible p = 1 is 
small, limiting the room for improvement. Alternatively, 
I c can be reduced by decreasing H e ff, which at small 
external field H is dominated by the anisotropy of F2. 
Simply reducing the magnetization M 2 of F2 would de- 
crease H e f f , but this would also compromise the stability 
of the magnetic configuration. Devices with perpendic- 
ular magnetic anisotropy can in principle overcome this 
limitation, but the potential advantages were offset by 
the reduced polarization and increased a Attempts 
to reduce m 2 while stabilizing F2 with an antiferromag- 
net have encountered similar complications [f|. 
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FIG. 1: Dynamics of two magnetic moments mi and m2 due 
to the spin torques rf TT and r^ TT exerted on the respective 
moments, (a) the signs of gi and <?2 are the same, (b) the 
signs of g\ and gi are opposite. 



Here, we discuss the previously unexplored mecha- 
nism affecting the efficiency of ST, involving simultane- 
ous current-induced effects on both ferromagnets in a bi- 
layer consisting of Fi and F2. This mechanism couples 
the dynamics of the two magnets. We describe analytic 
results for a simple model system, and present realistic 
numerical calculations in macrospin approximation. Our 
results demonstrate that I c can be either increased or 
reduced with respect to the value given in Eq. [T] by a 
suitable choice of Fi and F2. Most importantly, we show 
that the dynamics of both magnets are always excited si- 
multaneously, and thus both magnetic layers always par- 
ticipate in the current-induced behaviors. 

Analytic Model. To introduce the idea of dynamical 
coupling of magnetic layers by STT and understand the 
consequences of such coupling, we consider an idealized 
system of two nanoscale magnets Fi and F2 represented 
by the magnetic moments mi and m 2 , whose magnitudes 
are fixed in the macrospin approximation. For simplicity, 
this model neglects the demagnetizing fields of the fer- 
romagnets and their dipolar coupling. The dynamics of 
the magnetic moments in the presence of current I > 
flowing from Fi to F2 can be described by two Landau- 
Lifshitz equations coupled by the STT term Q 



dmi, 
dt 



= jmx t 2 x [H- 
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where 7 is the gyromagnetic ratio, and Si t 2 = ^1,2/(7^) 
are the spins of Fi,F2. For simplicity, we assume the 
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same a for both magnets, and neglect the dependence of 
gi.2 on the relative orientations of mi and mi. The exper- 
imentally determined values of I c are similar in both par- 
allel (P) and antiparallel (AP) configurations of the mag- 
netic moments, justifying the latter approximation [4|. In 
the limit Si 3> S2, the last term in Eq. [2] for mi is neg- 
ligible due to the large value of Si, resulting in a static 
solution for ni\. The equation for mi then yields preces- 
sion of m,2 at / = I c o determined by Eq. [TJ consistent 
with the models for the dynamics of a single ferromag- 
net. At / > I c q, the AP state becomes stable, and the 
solution for m-2 is a precessional reversal into this state. 

We are interested in the solution of Eq. [2] for com- 
parable values of Si and S2. As illustrated in Fig. [TJ 
STT exerted on both layers can induce their simultane- 
ous precession around the magnetic field H with cone 
angles 9\ and 62- A simple estimate shows that a sta- 
ble configuration involves precession of mi and mi with 
the same polar angle (Fig. [lja) ) . To demonstrate the 
stability of this configuration, we assume that mi lags 
behind m2 by a small polar angle </>. The torque t stt \ 
exerted on mi then acquires an additional component in 
the polar direction, resulting in the increase of its angu- 
lar frequency by Au>i = jIgisin(9i)sin(4>)/[eSisin(9i)]. 
Similarly, a component of STT exerted on mi in 
the polar direction increases its angular frequency by 
Aui = ^flgisin{6i)sin(4>) /[eSisin(Qi)]. The stability 
of the coupled precession with respect to fluctuations 
of <j) requires that Awi > A.U1, which is satisfied if 
[sin{6\) / sin{6i)] 2 < Si/ Si. The derivation given below 
shows that the latter inequality usually holds, ensuring 
the stability of in-phase precession. 

The stationary form of Eq. [2] for in-phase precession is 

Igisin(9i — 9\) = eS\aH sin(9i) 
Igisin(6i — 9\) = eSiaHsin(9i). (3) 

The filtering properties of Fi and F2 are often similar, 
resulting in gi rj gi. Eqs. [3] then lead to 

Sism(^) = Sisin(6 2 ) 7 (4) 

which satisfies the condition for the stability of in-phase 
precession, as discussed above. Eqs. [3] yield the current 
for the precession onset I c = 7 c o/[l — S2/S1]. We em- 
phasize three main consequences of Eqs. [3] Firstly, both 
moments always precess simultaneously. Secondly, the 
relation between the amplitudes of the dynamics of the 
two moments is determined by the ratio A = Si/ Si. In 
particular, the dynamics of mi is negligible at A <C 1. 
Precession of mi can then occur only at I = I c > 0. 
Conversely, the dynamics of mi is negligible at A ^> 1, 
while precession of mi can occur at / < 0. Finally, I c 
also depends on A. In particular, it diverges when the 
magnetic moments become equal (A=l). 

To determine the nature of excitations and stability 
regimes for the magnets coupled by STT, we insert Eq.0] 
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FIG. 2: Stability diagram for the analytic model of two-layer 
system with varied A — S1/S2. (a) gi = (72, bistable regions 
are hatched, (b) gi = —gi, the region of stable precession is 
filled. For A < 1, 7 c0 is defined by Eq. [TJ For A > 1, J c0 is 
defined by Eq. [JJwith rri2 replaced by mi to emphasize the 
dominance of this layer's dynamics. 

into the second of Eqs. [3j leading to 



^Jl/A 2 - 1 + cos 2 {9 2 ) - cos(6 2 ) ' 

where I c q — aHeSi / gi . Eq. [5] describes a monotonically 
decreasing function of 9, which implies that precession 
is unstable at any /. On the other hand, both 9i = 
and 62 = 7T are stable in the range I' c = j^j < I < 
= I c , where I' c is the current at which the AP state 
becomes stable. The resulting stability diagram is shown 
in Fig.[5Ja), where the bistable regions are hatched. 

As a consequence of the dynamical coupling due to 
STT, the value of I' c is reduced by a factor of two at 
A = 1, while I c diverges. It remains to be seen whether 
these coupling effects can be utilized to increase the ef- 
ficiency of magnetoelectronic devices. Current-induced 
bistability in devices with A « 1 should generally result 
in telegraph-type noise due to thermally activated transi- 
tions between the two stable configurations. Such noise is 
detrimental to most applications. However, such config- 
uration may be useful if small-amplitude current-induced 
dynamics is undesired. For A < 1, dynamical coupling 
also makes it impossible to induce dynamical states at 
/ < 0. Similarly, dynamics cannot be induced at / > 
for A > 1. We shall see below that this result also holds 
for a more realistic model. The fact that bipolar current- 
induced excitations were observed [1, 0] must indicate a 
significant breakdown of the macrospin approximation in 
these experiments. 

Applications of magnetoelectronic devices in mi- 
crowave generation require I c to be minimized, and sta- 
ble precession to be achieved over a significant range of 
/. We demonstrate below that I c can be reduced by 
the current-induced coupling when g\ and gi have op- 
posite signs. This requires reversing the spin anisotropy 
of one of the ferromagnets, which can be accomplished 
by doping Fi or F2 with impurities providing appropri- 
ate spin-dependent scattering in their bulk, and/or by 
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choosing N that inverts the spin anisotropy of electron 
scattering at F/N interface [10|. For simplicity, we now 
assume that gi = —g 2 < 0. The torques exerted on the 
two magnets now result in the mutual attraction of the 
moments at I < 0, resulting in a stable collinear configu- 
ration at 7 < for any A. Dynamical states are induced 
only by I > 0, regardless of the value of A. Fig. QJb) 
shows that the relative precession phases of mi and m,2 
are now shifted by 180°. Estimates similar to the ones 
presented above for gi, g 2 > show that this precessional 
configuration is stable. 

Eq.[4]for the precession cone angles of the two moments 
still holds, but Eq. [5] is replaced with 



Ida I A 



y/l/A 2 - l + cos 2 (0 2 ) + cos(9 2 ) 



(6) 



The resulting stability diagram shown in Fig. OJb) in- 
cludes a region of stable coupled precession of the two 
magnets. At A — 1, the precession onset current is re- 
duced by a factor of two, while the largest value of 7 at 
which the precession remains stable diverges. 

Numerical Simulations. To determine whether the 
current-induced coupling remains robust in realistic sys- 
tems, we solved Eqs. [2] numerically. We used the ge- 
ometry and the magnetic properties typical for the spin- 
transfer devices, in which Fi and F2 are Py layers with 
M = 800 emu/cm 3 and thicknesses t\ and i 2 , patterned 
into an elliptical shape with dimensions 120 x 60 nm 
We used the material parameters known from the magne- 

The demagnetizing fields 
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totransport measurements 
were taken into account. The dipolar coupling was ne- 
glected to eliminate its influence on the interpretation 
of current-induced behaviors. We included the depen- 
dencies of the spin transfer efficiencies gi,g2 on the con- 
figuration of the magnetic system to take into account 
the slight experimental asymmetry between the P^AP 
and AP— >P reversals The equations were numerically 
integrated by the stochastic Heun method with a fixed 
time step set to 2 psec. To verify the convergence, the 
step size was decreased by a factor of two, which did not 
significantly affect the results. Random field approxima- 
tion was used to model the thermal activation between 
different current-induced modes, with temperature set to 
50 K. Damping parameter a = 0.03 was used for both 
magnets [12j . 

The model was tested for three limiting cases: t\ 3> t 2 , 
ti <C t 2 , and ii = t 2 . In the first case, STT induced the 
dynamics of only m 2 , and only at J > 0. Fig. [3]^ a) shows 
the normalized spectra of the component of m 2 along 
the hard in-plane axis for A = 0.2. The peaks corre- 
spond to several harmonics of the precession frequency. 
After the onset of small-angle precession at I c — 1.2 mA, 
the frequency of precession decreases due to transition to 
large-angle clamshell precession trajectory. The transi- 
tion to the out-of-plane precession at 7 > 2.6 mA results 
in the increase of precession frequency. These results 
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FIG. 3: (a)-(c) Normalized spectral intensity for the dynamics 
of the in-plane hard-axis component of 7712, (a) A = 0.2, (b) 
A — 0.67, (c) A = 1 and gi — —gi. Lighter colors correspond 
to higher intensity. The same scale is used in all three plots, 
(d) Time trace for the trajectories of mi (dark blue) and ni2 
(light pink) for A = 0.67, I = 15 mA. x is the easy direction, 
y is the hard in-plane direction, and z is perpendicular to the 
film plane. The calculations were performed at H — 1 kOc 
along the easy axis, and a fixed t2 = 4 nm. 



are consistent with the published calculations for single 
layer dynamics, and are in overall agreement with the 
experiments |13j . Similar spectra were obtained for the 
dynamics of mi at / < 0, for t\ <C t 2 . Calculations for 
ti = t 2 yielded identical spectra for mi and m 2 regardless 
of the sign of /, as expected for this symmetric geome- 
try Additionally, the calculations for t\ — t 2 and H = 
reproduced the sequential thermally activated flipping of 
mi and m 2 experimentally observed in symmetric mag- 
netic nanopillars 

Current-induced coupling of the dynamics of two lay- 
ers significantly affected the spectra at 0.4 < A < 1, as 
illustrated in Fig. Efb) for A — 0.67. The onset current 
for the magnetic dynamics is larger than in Fig. [3Ja), 
in agreement with the results of the analytic model pre- 
sented above. At 3 mA< I < 6 mA, the data exhibit 
a broad incoherent excitation spectrum due to random 
transitions between different types of precession dynam- 
ics of m 2 . Specifically, the precession of m 2 alternates 
between the clamshell-type and the out-of-plane trajec- 
tories, depending on the relative phase with the elliptical 
precession of m 2 . The latter is not phase-coherent with 
7712. The formation of a sharp peak at 6 mA< / < 12 mA 
is associated with the complete transition to out-of plane 
precession of m2, resulting in nearly static deflection of 
mi towards m 2 . At 7 > 12 mA, the oscillations of mi 
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become sufficiently large to disrupt the periodic preces- 
sion of ra2, resulting in increasingly chaotic dynamics of 
both. Fig.[3^b) shows the trajectories of mi and m 2 over 
a 0.5 nsec period at 15 mA. Despite nearly chaotic be- 
havior of both moments, some parts of the trajectory of 
mi resemble clamshell and out-of-plane precession. The 
chaotic behaviors persisted in deterministic simulations 
at T = 0, eliminating thermal fluctuations as their cause. 
Based on the analysis of the trajectories, we believe that 
these behaviors are caused by a combination of large 
phase space associated with the dynamics of both mag- 
nets, and nonlinear nature of the large-amplitude mag- 
netic dynamics. 

At A > 0.67, the broad excitation band at small / and 
the out-of-plane precession peak gradually disappeared. 
At A = 1, only the I > 12 mA continuum remained, and 
similar features appeared at I < —12 mA. There were 
no excitations at / < for all A < 1, in agreement with 
the analytic results described above. These results seem 
to be inconsistent with the intuitive picture, according 
to which / > can induce the dynamics of m,2, while 
/ < can induce the dynamics of mi even for A < 1 [l5j |. 
However, analysis of the magnetic trajectories shows that 
all oscillations of mi at I < are efficiently suppressed 
by W2 closely following mi, thus reducing STT exerted 
on both layers. 

Calculations for g\ — —g% showed that I c is reduced in 
this case, in agreement with the analytic model. Fig. [3fc) 
shows spectra calculated for A = 1. I c = 0.75 mA is 
slightly over half of the value I c o obtained for A <C 1. 
Therefore, reduction of the critical current for gi = —g 2 
appears to be a robust feature of dynamical coupling. 
The spectra do not exhibit broadening and chaotic dy- 
namics characteristic of the g\ — g 2 data. The peaks in 
Fig. [3Jc) abruptly terminate at I = 3.8 mA due to the 
formation of a static stable configuration with mi and m 2 
oriented opposite to each other perpendicular to the film 
plane. Such static configuration is, however, not present 
for A < 1. 

We have performed additional calculations without 
some of the simplifying assumptions of our model. In par- 
ticular, we tested the effects of different magnetic damp- 
ing in the two magnets and different magnetizations. We 
also checked the effect of the dipolar coupling between the 
layers. The modifications did not qualitatively change 
our conclusions regarding the effects of current-induced 
coupling, as long as the resonant frequencies of the two 
layers remained similar to each other. 

In summary, we have analyzed the simultaneous effects 
of spin torque on both layers in a magnetic bilayer. We 
showed that the spin torque results in dynamical coupling 
between the layers, modifying their individual current- 



induced dynamical properties. In particular, the onset 
of magnetic precession of both magnetic always occurs 
at the same current. In case of similar spin-transport 
anisotropy of both layers, the onset current for magnetic 
precession is increased, diverging for identical values of 
the moments of the two layers. Realistic numerical sim- 
ulations show that coupling often leads to incoherent dy- 
namical regimes associated with transitions between dif- 
ferent dynamical states of the bilayer. However, for the 
opposite spin-transport anisotropy of the two layers, the 
onset current for the magnetic precession is decreased, 
and coherence is maintained even for the identical di- 
mensions of the magnetic layers. This effect of dynamical 
coupling on the spin transfer efficiency can become useful 
for implementing devices with improved performance. 

This work was supported by the NSF Grant DMR- 
0747609. I thank Lidia Novozhilova for help with numer- 
ical simulations. 
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